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ABSTRACT 

We show that in the framework of the pinch technique the univer¬ 
sal part of the p parameter can be meaningfully defined, beyond one 
loop. The universal part so obtained satisfies the crucial requirements 
of gauge-independence, finiteness, and process-independence, even when 
subleading contributions of the top quark are included. The mechanism 
which enforces the aforementioned properties is explained in detail, and 
several subtle field theoretical issues are discussed. Explicit calculations 
of the sub-leading two-loop corrections of order 0(G 2 M^) are car¬ 

ried out in the context of an SU(2) model, with Mw = Mz } and various 
intermediate and final results are reported. 
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1. Introduction 


One of the most important quantities of the standard model is the p-parameter [1], de¬ 
fined as the ratio of the relative strength between neutral and charged current interactions, 
at low momentum transfer, namely 

_ Gnc(O) _ 1 ( , 

Gcc( 0) 1 — A p 

The p-parameter displays a strong dependence on m t and affects most electroweak param¬ 
eters such as Ar, M\y , and sin 2 9 e ff(Mz)- The p parameter defined above as the ratio of 
two amplitudes is a gauge independent and finite quantity. In addition, it is manifestly 
process dependent, and its value depends on the quantum numbers of the external parti¬ 
cles chosen. To fully determine the value of p for a given neutral and charged process, one 
must compute the complete set of Feynman diagrams ( self-energy, vertex and box graphs) 
to a given order in perturbation theory. However, traditionally one focuses instead on the 
quantity A, defined as 


_ A w ( 0) _ A z ( 0) + (2s/c)A lZ (0) , , 

M* w Ml ’ C ' 

where Aw, A z and A~. f z are the cofactors of g^ in the WW and ZZ self-energies, respec¬ 
tively, and s = sin 6w and c = cosOw- The parameter A is often called the “universal” 
part of p, since, by definition, does not depend on the details of the process. According to 
the standard lore A contains the dominant contribution to p. 

In this paper we address theoretical issues related to the definition and calculation of 
the universal part of the p parameter. Several of the ingredients of the subsequent analysis 
are known, but they exist in a fragmented fashion and their relevance in the context of p- 
parameter calculations has not been fully recognized. The purpose of the present work is to 
provide a unified view of all relevant facts, incorporating them into a coherent framework. 
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The leading m t contributions to A, both from one-loop and two-loop diagrams, are 
gauge-independent and ultraviolet finite However, as soon as the subleading contribu¬ 
tions are considered, A becomes gauge-dependent and ultraviolet divergent. More specif- 
icaly, the leading one-loop m t contributions to A (of order G^mf) are manifestly gauge- 
independent, since the gauge fixing parameter does not appear inside fermion loops, and 
ultraviolet (UV) finite. On the other hand, as was explicitly shown by Degrassi, Kniehl, 
and Sirlin [2], the one-loop bosonic contributions to A (subleading in m^, of order j 2 mJ) 
are gauge dependent and, except when formulated within a restricted class of gauges, UV 
divergent [3]. Similarly, the leading two-loop contributions to A [4-6](of order G^rnf) 
are also gauge-independent and UV-finite, exactly as their one-loop counterparts. On the 
other hand, the straighforward calculation of the subleading two-loop m t contributions (of 
order G'^m^M'^) to A, carried out in the renormalizable Feynman gauge = iz = 1), 
gives rise to an answer which is ultraviolet divergent, in the sence that y terms survive 
[7-8]. In addition, as stated in [7], the result for A is gauge dependent, in the context of 
the R^ gauges. The reason for this rather striking analogy between the one-loop and the 
two-loop analysis, is the fact that, as soon as the subleading contributions (of order mj 
at one-loop, and of order at two-loops) are considered, Feynmam graphs containing 
the W and Z gauge-bosons must be included. It is the inclusion of such graphs, which, 
when carried out without a concrete guiding principle, gives rise to the aforementioned 
pathologies. 


In order to understand the origin of the problems associated with the subleading 
contributions one has to first establish the mechanism which enforces the good behavior 
of the leading contributions, in particular their UV finiteness. If we denote the leading 
contributions (both at one and two loops) to the to WW and ZZ self-energies by [n^]^ 
and [n^]^, respectively (the superscript l stands for ’’leading”), and use the fact that 
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A^z(O) = 0, we can write for A^: 


^(7) _ _ ^-^(0) 


It is elementary to show (see section 3) that is finite as long as the following relation 


holds 


yW yW 

L W L Z 


i_ fsf\ (e) 

/ /7\ I l O I 


where Zw and Zz are the wave function renormalization constants of the W and Z fields, 


& 2 V 0 _ Re[A { ‘\MD) 


respectively, and 


Eq. (1.4) is indeed true for the leading contributions, due to a set of QED-like Ward 
identities relating the vertex and fermion wave function renormalization constants (the 
QED analog of Z\ = Z 2 ). An equivalent way for establishing the finiteness of the leading 
contributions is to resort to the following Ward identities: 


q»qA n 


/^iG) _ a/T 2 ttW 


= 


q,qA^zf ] = , 

where 11^, and II X are the <(><(> and yy self-energies, respectively, with ^ (y) the charged 
(neutral) would-be Goldstone bosons. Writing the general self-energy in the form 


U-tiv^q ) — A(q + B(q ) 




and using the algebraic identity 


• 4 ( 0 )= -A{q“q" n,.„(<i 2 )} . 


together with the Ward identities of Eq. (1.6), and the fact that A^(0) = 0, we can write 
for A^: 

^(7) _ -^ z \ 0) 




4 



The final ingredient which enforces the finiteness of the result for the leading A p con¬ 
tributions, is the equality = Z^\ where Z^ and Z^p are the wave-function renor¬ 
malization constants of and Hx\ respectively (see section III). However, none of the 
above conditions are valid anymore, when one calculates the subleading parts of the W 
and Z self-energies in the framework of the R £ gauges. Consequently, since the mechanism 
enforcing the finiteness does not operate any more, the resulting expressions do not have 
to be UV finite, and, indeed, they are not. 

The standard way to circumvent the above problems is to readily abandon the notion 
of a “universal” part of Ap, by stating that, instead of only A, the entire process must 
now be considered, in order to restore the finiteness and gauge-independence of the final 
answer. So, one has to introduce vertex and box corrections, which render the result 
gauge-idependent and finite, at the expense of making it process-dependent, and therefore 
non- “universal”. 

This unpleasant trade-off between gauge-independence and process-independence can 
be avoided however [2] if one defines the universal part of p within the framework of 
the pinch technique (PT) [9-15]. As shown in Ref [2], the PT at one-loop gives rise to 
a gauge-independent and UV finite answer, vjithout introducing any process-dependence. 
This is so because all PT self-energies are individually gauge-independent, and in addi¬ 
tion, all conditions which enforce the finiteness are valid for both leading and subleading 
contributions. Evidently, the PT restores the mechanism for the cancellation of the UV 
divergences, and at the same time guarantees the gauge-and process-independence of the 
final answer. 

In this paper we propose to elevate the A defined in terms of the PT self-energies as 
the truly universal part of Ap, beyond one-loop, in the spirit originally suggested in Ref [2]. 
This new quantity, which we denote by A, is endowed with three crucial properties: 

i) A is independent of the gauge-fixing procedure and the gauge-fixing parameter, 

ii) A is ultraviolet finite, and 
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iii) A is process-independent. 

In addition to the above obvious theoretical advantages, the calculation of A is sig¬ 
nificantly fascilitated by the fact that only self-energy-like graphs need be considered. In 
particular, no vertex or box diagrams need be calculated, in order to render the answer 
finite, as is the case in the conventional treatment [16]. Therefore A lends itself as the 
natural generalization of the conventional A, which can consistently accomodate both lead¬ 
ing and subleading contributions, and can be expressed in a closed analytic form up to two 
loops. 

In Ref [7] the nature of the two-loop subleading contributions, as well as their nu¬ 
merical importance for the conventionally defined A, were studied in the context of an 
SU( 2) model, with Mw = Mz = M, and no photon. In Ref [8] the previous analysis 
was extended to the full standard model; the two answers turned out to be numerically 
rather close. Since we are mainly interested in addressing the conceptual issues involved, in 
this paper we also restrict ourselves to the study of this simplified version of the standard 
model. Of course, we have no a priori knowledge if the A of the SU( 2) and the A of the 
full stadard model will be numerically close, as was the case between the results of Ref [7] 
and Ref [8]. 

The paper is organized as follows: In section 2 we briefly review the PT, mainly as it 
applies in the present context. In section 3 we define the universal quantity A and discuss 
some of its properties. In section 4 we calculate the subleading top quark contributions to A 
in the framework of the SU(2) model. This section is rather technical and contains several 
intermediate results. Finally, we present our conclusions in section 5. In addition, we 
present two Appendixes concerning scalar two loop integrals at zero momentum transfer, 
and the Feynman rules we have used. 
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2. The pinch technique 


Gauge-invariant effective self-energies 

The simplest example that demonstrates how the PT works is the gluon two point 
function. [10] Consider the S'-matrix element T for the elastic scattering such as qiq 2 —» 
qiq 2 , where qi,q 2 are two on-shell test quarks with masses m\ and m 2 . To any order in 
perturbation theory T is independent of the gauge fixing parameter £. On the other hand, 
as an explicit calculation shows, the conventionally defined proper self-energy depends on 
£ defined through the tree level gluon propagator 

= • ( 2 - 1 ) 

At the one loop level this dependence is canceled by contributions from other graphs, 
which, at first glance, do not seem to be propagator-like. That this cancellation must 
occur and can be employed to define a g.i. self-energy, is evident from the decomposition: 

T(s,t,m 1 ,m 2 ) = T 0 (f,£) + Ti(f,m l5 £) + T 2 (f,m 2 ,£) + T 3 (s, t, m 1 , m 2 , £) , ( 2 . 2 ) 

where the function To(t,£) depends kinematically only on the Mandelstam variable t = 
— (pi — pi ) 2 = ~q 2 , and not on s = (p\ ~\~ P2) 2 or on the external masses. Typically, 
self-energy, vertex, and box diagrams contribute to To, Ti, T 2 , and T 3 , respectively. Such 
contributions are £ dependent, in general. However, as the sum T(s, t, mi, m 2 ) is g.i., it is 
easy to show that Eq. (2.2) can be recast in the form 


T(s,t,mi,m 2 ) = T 0 (t) + Ti(t,mi) + T 2 (t,m 2 ) + T 3 (s,t,mi,m 2 ) , (2.3) 

where the Ti (i = 0,1, 2, 3) are individually ^-independent. The propagator-like parts of 
vertex and box graphs which enforce the gauge independence of To(t), are called pinch 
parts. They emerge every time a gluon propagator or an elementary three-gluon vertex 
contributes a longitudinal k^ to the original graph’s numerator. The action of such a term 
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is to trigger an elementary Ward identity of the form ft = (j) + ft — m) — (j) — m) when it 
gets contracted with a 7 matrix. The first term removes (pinches out) the internal fermion 
propagator, whereas the second vanishes on shell. From the g.i. functions Ti (i = 1,2,3) 
one may now extract a g.i. effective gluon (G) self-energy fl^^lq), g.i. Gqiqi vertices T^\ 
and a g.i. box B , in the following way: 


T 0 = g 2 u 1 j^u 1 [(\)fl tlu (q)(\)\u2j l 'u 2 , 
q z q z 

Ti = g 2 u 1 (\)u 2 Yu 2 , 

r (2.4) 

T 2 = g 2 uij^u 1 (^)u 2 f'i 2) u2 , 
q z 

Ts=B , 

where Ui are the external spinors, and g is the gauge coupling. 

The one-loop expression for calculated in the Feynman gauge £ = 1 is given 

by [ 10 ]: 


n F ( g ) = n« =1 )( g ) + vn p ( g ) , ( 2 . 5 ) 

where t^y = —g^y + q^q 1 '/q 2 , l^y = q^q 1 ' /q 2 are the usual transverse and longitudinal 
projectors, and 

nF(q) = - tc ^ q2 S D nk + q r- ■ (2 ' 6) 

where j D = g 4 ~ D J * s the dimensionally regularized loop integral, D is the dimen¬ 

sionality of space-time, and c a is the quadratic Casimir operator for the adjoint represen¬ 
tation [for SU(N) } c a = N], After integration and renormalization we find 

n p ( 9 ) = 2 c „ (l |b7ln(^T)] . (2.7) 


Adding this to the Feynman-gauge proper self-energy 



we obtain for II ^y(q) 

-q 2 

hl/tj/^) — bg ln( — )tnv , 

g 

8 


( 2 . 8 ) 


(2.9) 



where b = 16 1 7r 2 is the coefficient of — g 3 in the usual /3 function of QCD without 
fermions. This procedure can be extended to an arbitrary n-point function; of particular 
physical interest are the g.i. three and four point functions [17-18] and T^ap [19]. 

Finally, the generalization of the PT to the case of non-conserved external currents is 
technically more involved, but conceptually straightforward [20-21], 

The current algebra formulation of the pinch technique 


An important alternative formulation of the PT in the context of the SM has been 
introduced by Degrassi and Sirlin Ref. [12]. In this approach the interaction of gauge bosons 
with external fermions is expressed in terms of current correlation functions, i.e. matrix 
elements of Fourier transforms of time-ordered products of current operators [22], This 
is particularly economical because these amplitudes automatically include several closely 
related Feynman diagrams. When one of the current operators is contracted with its four- 
momentum (i.e. the four momentum absorbed by the current), a Ward identity is triggered. 
The pinch part is then identified with the contributions involving equal-time commutators 
in the Ward identities, and therefore involve amplitudes in which the number of current 
operators has been decreased by one or more. As emphasized in Ref. [12], this procedure 
has an important advantage when one considers external particles endowed with strong 
interactions. Because the contributions from the equal-time commutators are not affected 
by the dynamics of the strong interactions, the aforementioned identification ensures the 
universality of the ’’pinch parts”. That is, the cofactors of the current operators in the 
pinch parts are the same whether the external particles are leptons or strongly interacting 
fermions. To illustrate the method with an example, consider the vertex function iU^(W) 
that contributes to Fig.(lb), where now the gauge particles in the loop are W’s, the 
incoming one is a Z, and the incoming and outgoing fermions are massless. It can be 


9 



written as 


iUhW) 


ig 3 c 


d n k 


A'lik) A'L(k + q)W‘°( 2k + ,) A - 9 A "(2 q + kf- g“\k - q)’} 


2 J (2x)" '"’ tf 

x f d'kxe*’ < f\T'\J$(x)J$ v (Q)\\i >■ 


( 2 . 10 ) 


where 

(2 ' n) 

with i = W,Z, 7 and M~ i = 0, are the propagators of the gauge bosons in a general R^ 
gauge. An appropriate momentum, say k a , from the three gauge boson vertex or the 
longitudinal part of the propagator can be transformed into a derivative ^— acting on 
the T* product. Invoking current conservation this leads to an equal-time commutator of 
current operators. Thus, such contribution are proportional to the matrix element of a 
single current operator, namely < >', these are precisely the pinch parts. Calling 

iU^(W) p the total pinch contribution from Eq. (2.10), we find in the £ = 1 gauge 


U*(W) p = ig 3 c<f\J*\i> 


i 

[V - Mi)[[k + qf - Ml] 


( 2 . 12 ) 


Clearly, the integral in Eq. (2.12) is the generalization of the QCD expression Eq. (2.6) to 
the massive gauge boson case. 

Ward identities of the PT 

Another important fact is that the PT Green’s functions satisfy tree-level Ward iden¬ 
tities. Most noticeably, the g.i. QCD vertex Ggigi satisfies the following Ward-identity: 


g^Tfq = E(p + q) - £(p) (2.13) 

where E(p) is the g.i. quark self-energy [23] The above QED-like Ward identity, which 
is not true for the conventional enforces the equality Z\ = Z 2 between the vertex 

renormalization constant Z\ and the quark wave function renormalization constant Z 2 . 
Consequently, exactly as happens in QED, the PT vacuum polarization contains the entire 


10 



running of the QCD coupling, as shown already by the explicit result of Eq. (2.9) [10], 

[24]. 

The above QCD results have been generalized for the electroweak part of the standard 
model, where it was found that the one-loop PT Green’s functions satisfy again tree-level 
Ward identities [18], [21], [25]. Therefore, the wave function renormalizations for the 

PT 77 and WW self-energy contain the running of the gauge couplings e 2 (q 2 ) and g 2 (q 2 ), 
respectively [12], [13]. Denoting by fl^ the gauge boson PT self energies (f[^ or ft^), 
by 0^ the mixed PT self energy of a gauge boson and its associated unphysical scalar 
(nJT* = n+ = -n;. or and by O the PT self energies of the unphysical scalars 

(fl^ or fl x ), the following WI hold [26]: 

q^Tl^ - iMQ„ = 0 

q II Q^ + iMQ = 0 (2.14) 

q’ M q’'% v — M 2 0 = 0 . 

Additional WI between other PT Green’s functions can be found in the literature. As 
was explained in detail in [21], the PT Ward identities are instrumental for the final 
cancellation of gauge dependences in S-matrix elements. 

Imposing the elementary requirement that the renormalized PT Green’s functions 
should respect the same Ward identities as their unrenormalized counterparts, we obtain 
the following relationships between the standard model renormalization constants: 


and 


with 


Zw = Z 


-2 


- - be 2 

Zt = 5 


A A A A , ^ M w 

ZjH — Zj x — Zjs — Lw -|-T— 

M 2 

m w 


8&_ _ SM (y 8M 2 Z 


(2.15) 

(2.16) 


(2.17) 


(2.18) 
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In the simplified SU(2) model we will use laetr we have: 


Zz = Zw = Z, 


-2 


(2.19) 


and 


Zh — Z x — Z<f, — Z a 2 + 


SM 2 
9 + M 2 


( 2 . 20 ) 


As usual, gauge boson self energies are cast in the form 


fi ^(q) = g^A(q 2 ) + ^B(q 2 ) 


( 2 . 21 ) 


and the renormalization constants are defined from the expansion 


A(q 2 ) = A(M 2 ) + (q 2 - M 2 ) dA ^ 


dq 2 


+ A(q 2 ) 


2 \ finite 


q 2 = M 2 


( 2 . 22 ) 


as 


SM 2 = Re ( A(M ) 


z- 1 = 1 - 


dA(q 2 ) 


dq 2 


q 2 = M 2 


(2.23) 


3. The universal A 


In this section we focus on the universal part of p defined by the PT. In particular we 
will emphasize issues of gauge-independence and finiteness. 

Traditionally the universal part A is defined as in Eq. (1.2). A vanishes in the limit of 
exact SU(2)v custodial symmetry, e.g. for Mw = AIz (no hypercharge) and for degenerate 
fermion doublets, m u = m& [27]. 

The fermionic one-loop contribution is given by [1] 


A p) = N c 


8a/2 tt 2 


m l + m 2 d + 


2m l m d 

m 2 — m 2 d 


i m d 

In—p 


(3.1) 
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Clearly, A p 1 -> 0 as m u —> mg- When the mass splitting in the fermion doublet is large, 
as in the case of the top and bottom quarks, the factor in square brackets in Eq. (3.1) 
is replaced by m 2 , the heavy fermion mass. By neglecting the contribution of all light 
fermions the one loop fermionic contribution is written as 


= N c x t 


(3.2) 


where 


G ,™ 2 


9 2 

g m t 


_ _ y fo o\ 

*' ~ 8V2^-~ 16x2 4M 2, (•>'•>( 

If one attempts to use the definition of Eq. (1.2) to include bosonic one loop corrections 
one is faced with two problems 

(a) The result is £ — dependent 

(b) Unless computed in a special class of gauges, it is ultraviolet divergent. 

In particular, regarding the first point, the dependence on the gauge fixing parameter 
enters through the tree-level propagators for the IU, the Z and the photon which, in the 
R g gauges they are given by Eq. (2.11). In addition, the tree-level propagators of the 
unphysical Goldstone bosons are given by 




(3.4) 


f/2 - f,-M? ’ 

with (s } i) = (<(>, IU) or (X,Z), and they also explicitly depend on The conventional 
one-loop self-energies depend explicitly on the gauge-fixing parameters, even at q 2 =0. 
Regarding point (b), unless the relation 


iw = i £jsin 2 9 + ^zcos 2 6 . 


(3.5) 


between the gauge fixing parameters ^ is satisfied, the resulting expression for A contains 
a term proportional to U 

The problems mentioned above persist when one computes the two-loop contributions 
to A. Again, as happens in the one-loop case, the leading contributions are both gauge- 
independent and finite. As soon as the subleading contributions are taken into account the 
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pathologies familiar from the one-loop, reappear: The results are again gauge-dependent 
and, even when computed in the Feynman gauges f^w = tz = C 7 = which obviously 
satisfies the one-loop condition of Eq. (3.5), are ultraviolet divergent. Evidently, Eq. (3.5) 
breaks down beyond one loop. 

Before we proceed to the study of A defined via the PT, it is worthwhile to fur¬ 
ther elaborate on the gauge independence of the leading two-loop contributions to the 
conventional A mentioned above. The reason for it is simply that the S-matrix is gauge- 
independent, and there are no vertex or box contributions proportional to which could 
cancell any possible gauge-dependences coming from the self-energy graphs. Consequently, 
one is allowed to choose any convenient gauge for calculating these leading self-energy con¬ 
tributions. In particular, in the Feynman gauge (£ = l)only graphs with scalars and 
fermions contribute to this order. This is not generally true however for an arbitrary value 
of the gauge fixing parameter £. The graphs of Fig.9 for example will give rise to leading 
mf contributions, due to the longitudinal parts of the gauge boson propagators. Obvi¬ 
ously the characterization of individual Feynman graphs as ” leading” and ” subleading” 
is a gauge dependent statement. When all relevant graphs are correctly accounted for, 
they will indeed conspire to furnish a unique ^-independent answer, which will clearly be 
identical to the one obtained in the Feynman gauge. It is instructive to briefly highlight 
the mechanism enforcing the cancellations of the gauge-dependences. To that end we can 
employ the elementary algebraic identity 

1 1 (1-QM 2 

q 2 -iM 2 q 2 - M 2 + (q 2 - M 2 )(q 2 - £M 2 ) 

in the graphs containing scalars and fermions. The first term in r.h.s. of Eq. (3.6) is the 

Feynam gauge scalar propagator, whereas the second term resembles the longitudinal part 

of the corresponding gauge boson propagator. The final cancellation proceeds after using 

the elementary Ward identity 

k^P L = PPl = S-^p + k)P L - PrSJ\p) + m t P L - mj P R , (3.7) 


3.6 
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where Pr^l = (1 ± 75 )/ 2 , triggered by the longitudinal term k^k 1 ' of the gauge boson 
propagator. It is this Ward identity, which when applied for both AT and k u , gives rise 
to an extra power of m 2 , thus converting pieces of a diagram, which is subleading in the 
Feynman gauge, into leading. 

If one defines A instead in terms of the effective WW and ZZ propagators obtained 
via the PT all problems associated with the gauge-independence and finiteness of the 
subleading parts are automatically solved. We denote A the universal part of A p defined 
via the PT as follows 


A = 


Aw(0) Az{ 0) 


M?v 


Ml 


(3.8) 


It is important to notice the absence of the jZ mixing term in the above definition; this 
is so because in the PT the jZ self-energy vanishes at q 2 = 0, e.g. fl^(0) = 0. The 
PT self-energies are individually independent of the gauge-fixing parameters, and when 
combined according to Eq. (3.8) they give a UV finite answer. 

Although the gauge-invariance of the result in the context of the PT is guaranteed 
by construction, its finiteness may be less obvious. There are two equivalent ways of 
understanding why the PT definition gives rise to a finite expression, both relying on the 
Ward identities presented in section 2. Writing the WW and ZZ self-energies in the form 


and 


Eq. (3.8) yields 


Aw(q 2 ) — Aw(M^y) + ( q 2 — M$y)[ 1 — ^w \ + Ayy(q 2 ) 


Az(q 2 ) — Az{M 2 z ) + (q 2 — Mf)[ 1 — Z z x \ + A z (q 2 ) , 


A A \d t v + M 2 


A vv(0) A z (0) 


w 


Ml 


(3.9) 


(3.10) 


(3.11) 


A , which contains the terms proportional to ^ ( and possibly finite pieces, which we 


neglect at this point), is given by 

■lui.l/f, i Az(M 2 z ) 


^ div 


Mh 


M| 


+ (%w Z z ) 


(3.12) 


= 0 
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where in the last step we used Eq. (2.16). So, A of Eq. (3.11) is finite for leading, sublead¬ 
ing, and bosonic contributions. 

Another way to establish the finiteness of A is the following: After computing the 
off-shell WW and ZZ self-energies, which have the form of Eq. (1.7), we use Eq. (1.8) 
together with the Ward identities of Eq. (1.6), we can write for A: 

i _Aw(0 ) Az(0 ) 

~~K, W 

=[A {rl ^nq ]|^ 0 

We than use the fact that 

n„(, 2 ) = n 4 ,(m 2 v ) + tg 2 - at 2 ,)[i - z p] + ig(, 2 ) 

and 

n x (9 2 ) = n x (M|) + (q 2 - m|)[i - z- 1 } + n{(, 2 ) 

where 1 — Z~ 2 = ^\:' 2 1 , r = end 1 — Z ~ 1 = ^^/'-J' ir Taking the difference of 

the two scalar self-energies, and using the fact that = Z x = Z we have 

n*( 9 2 )-n,( 5 2 ) = [fi^(M^)-n x (Ml)]+(i-z _ 1 )[Ml-M^]+[fiJ( s 2 )-n{( s 2 )] . 0.14) 

The first two terms on r.h.s. of the last equation are proportional to but they are 
constant, independent of q 2 . Therefore, upon differentiation with respect to q 2 they vanish. 
The third term is g 2 -dependent and finite, and after diferentiating it with respect to q 2 
and subsequently setting q 2 = 0 we obtain the UV finite expression for A. Clearly, the 
above proof of finiteness does not depend on the choice of the renormalization point; so 
instead of expanding around q 2 = M^y and q 2 = M%, we can equally well expand around 
q 2 = / 1 2 and q 2 = , respectively. 

It is important to emphasize that all properties of the PT self-energies stemming from 
the PT Ward identities hold for the corresponding conventional self-energies computed in 


3.13 
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the background field method (BFM) [28], for every value of the gauge-fixing parameter 
fqf and £q, used for the quantum fields, and to all orders in perturbation theory [29]. 
Consequently, in the BFM the finiteness of A is true for every value of the gauge fixing 
parameter and to all loops [30]. The final answer is however not gauge-invariant. This 
is so because in the BFM the gauge boson self-energies depend in general on the choice 
of gauge-fixing parameters [23](these gauge dependent terms are however UV finite); this 
remaining gauge-dependence does not cancel when the difference of the WW and ZZ self¬ 
energies is formed, in order to construct A, already at the one-loop level. So, the one loop 
bosonic contribution defined via the BFM off shell <fxf> and yy, which are both 

explicitly ^Q-dependent, has the form 


^b BFM Zq=1 ^(Cq)- 


(3.15) 


I 5q =i = A * is the gauge independent bosonic PT result, given by 


A l = 


9 c 2 — 8c 4 + h(5c 2 — 6) o 

i +-- L ^ 2 )- 


c 2 (c 2 — h)( 1 — h) 


log(fc) , (3.16) 


where h = Mfj/m 2 , and Mu is the mass of the Higgs scalar. A 1(£q) carries explicitly the 
gauge parameter dependence (one sets for simplicity = fqf = £q). A 1(£q) vanishes at 
= 1, and when Mw = Mz (-s = 0), but is non-zero otherwise; its explicit expression 
has been reported in [31]. As one can see from Fig. 1, the gauge dependent term A K^q) is 
unbounded from above and below, and is numerically significant. We observe that although 
the BFM endows the Green’s function with the desirable theoretical properties, it fails to 
address the crucial issue of gauge-fixing parameter independence, as any other gauge fixing 
procedure for that matter. Nevertheless, it provides a convenient starting point for the 
implementation of the PT [23]. 
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4. Subleading top contributions to A in the SU( 2) model 


In this section we apply the formalism developed thus far to the case of an SU(2) 
model, which corresponds to the standard SU(2) X 17(1), with electromagnetism turned 
off. This means that s = 0 and c = 1, or equivalently Mw = Mz = M, and there is no 
photon. In such a model, a non-vanishing A comes only from the mass-splitting within 
a fermion iso-doublet; in particular, there are no genuinely bosonic contributions to A, 
since Mw = Mz- We ignore the contributions of all light fermions, and concentrate on 
the top quark contribution. Since at one loop A contains only fermionic contributions it is 
automatically gauge independent and UV finite. Obviously, at one loop the PT definition 
coincides with the conventional one. 

At two loops there are two kinds of contributions: 

(i) The leading, of order mj, which originate from graphs containing fermions and 
only scalars, without gauge bosons. 

(ii) The sub-leading, of order m^, which originate from the graphs of (i), if scalars are 
replaced by gauge bosons (W or Z). 

The leading contributions of a very heavy top quark to the conventionally defined A 
have been first computed in [4], in the limit where Mw = Mz = Mjj = 0 . The case of an 
arbitrary Higgs mass Mjj, but still Mw = Mz = 0, was computed in [5], [6], and the case 
Mw = Mz = M ^ 0 was presented in [7]. In the above calculations the Feynman gauge 
was used; as already explained in section 3, this convenient choice of gauge is legitimate, 
since the result is guaranteed to be gauge-independent and UV finite. Clearly, the PT and 
the conventional definitions are identical for the leading part of the calculation. 

The sub-leading top contributions to the conventionally defined A were first addressed 
in [7]; it was explicitly shown that the resulting expressions contain left-over terms pro¬ 
portional to 7 In addition, it was pointed out that this result, calculated in the renor- 
malizable Feynman gauge, was in fact gauge-dependent, and was correctly argued that 
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due to this theoretical shortcommings the inclusion of subleading corrections deprives the 
conventionally defined A of any physical meaning. In order to restore gauge invariance 
and finiteness several contributions from vertex and box diagrams were included; however, 
since no guiding principle such as the PT was followed, these contributions rendered the 
answer process-dependent. Furthermore, it was sugested that, since the sub-leading con¬ 
tributions cannot be defined in a process-independent way, the possibility of resumming 
them [32] should probably be abadoned. Since no closed expressions for two-loop vertex 
and box graphs exist yet in the literature, the process-dependent parts were calculated ap¬ 
proximately, up to order 0(M 2 / m 2 ), for the case of v^e scattering. The final conclussion 
was that the part of Sp that one extracts with their method, for the case of v^e scattering 
receives sizable corrections due to sub-leading top contributions. 

In the context of the PT all aforementioned pathologies are automatically bypassed. 
The answer is gauge-independent by construction, UV finite, and manifestly process- 
independent. We conclude therefore that in the context of the PT there is no limitation 
whatsoever in defining the subleading top contributions to the universal A. In particular, 
the necessary condition for attempting the resummation of the universal part of p, i.e. 
the process-idependence, is still valid. It turns out that the relative size of the subleading 
contributions compared to the leading is in accordance to what one would naively expect 
from a power series whose expansion parameter r = M 2 /m 2 is of the order of 1/4. 

We now proceed to the more technical aspects of our calculations. 

(a) It has been known for years that when computing the PT Green’s functions any 
convenient gauge may be chosen, as long as one properly accounts for the pinch contri¬ 
butions within that gauge [10]. In the context of the “renormalizable” R^ gauges the 
most convenient gauge-fixing choice is the Feynman gauge (£ = 1). This is so because 
the longitudinal parts of the gauge boson propagators, which can pinch, vanish for £ = 1, 
and the only possibility for pinching stems from the tree-boson vertices. As was recently 
realized [29], the task of the PT rearrangment of the S-matrix can be further fascilitated, 
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if one quantizes the theory in the context of the BFM. Even though the Feynman rules 
obtained via the BFM are rather involved, they become particularly suited for one-loop 
pinching, if one chooses the Feynman gauge (£q = 1) inside the quantum loops. In fact, all 
possible one-loop pinch contribution are zero in this gauge; consequently, the one-loop PT 
Green’s functions (which one can obtain for every gauge) are identical to the conventional 
Green’s functions, calculated in the Feynman gauge of the BFM. This property of the 
Feynman gauge in the BFM persists in two loop calculations, only for the subset of dia¬ 
grams which contain at least one fermion-loop [33], which is precisely the type of graphs 
we are interested in [34], We therefore choose to work in the BFM Feynman gauge. This 
correspondence between PT and BFM at = 1 [35] breaks down for the two-loop purely 
bosonic part [36]. The technical details leading to this conclusion will be presented in 
[37], 

(b) Using the algebraic identity of Eq. (1.8), and the WI of Eq. (2.14), we write A as 
in Eq. (3.13). Therefore, the entire calculation boils down to calculating the derivative of 
each Feynman graph shown in Fig.2-Fig.7, at q 2 = 0. In these figures we show the complete 
set of two-loop irreducible graphs that contribute only mf terms to A, in the SU( 2 ) model. 
Graphs with photons or with couplings proportional to sin 6w in the standard model are 
omitted. For example, graphs such as those shown in Fig .8 for the self-energies of the 
scalars do not contribute, since they contain couplings proportional to sin 6w = 0. In 
our calculation we have used a fully anticommuting 75 since this does not produce any 
inconsistencies for the graphs we have to compute. 

The validity of the Ward identity of Eq. (2.14) for the full standard model has been 
verified by directly contracting individual graphs by q^q 11 , before carrying out any loop 
integrations. In fact, the Ward identities of Eq. (2.14) hold individually for each of the 
graphs shown in Fig.2- Fig.7, where the corresponding graph of the gauge boson self energy 
-which is to be contracted by - can be obtained by replacing the external iff ( XX) 

legs by WW ( ZZ), respectively. The only exceptions are some of the triangle graphs, i.e. 
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Fig.5(a,c), that need to be combined with graphs such as those depicted in Fig. 8 , in order 
to yield zero in the r.h.s. of the WI. For example we have 


qVq 1 ' (5 ' a) + (8 ' C) ] - M w [ n ^ (5 ' a) + (8 ' d) } = 0 (4.1) 

\n^ (5 - c) + (8 - a) j - M^y \n^ (5 - c) + (8 - 6) j = o ( 4 . 2 ) 

In the above equations, we note that although the graphs of Fig .8 for the self energies of 
the scalars vanish for sin 6w = 0 , the corresponding ones for the gauge bosons do not. 
Using the notation 

Gk(q 2 ) = — i (Feynman — Graph k) 


we find it more convenient to act with the four-Laplacian instead of the regular derivative 


of q 2 . Namely 


dG k (g 2 ) 

dq 2 


2 0 “ 2D \P* Gk ( q 

q z = 0 \ / q =0 


This facilitates the computation enormously since it reduces it to straightforward alge¬ 


bra that can be carried out easily by hand. This procedure reduces each graph down to 


standard scalar two loop integrals at zero external momentum, for which explicit expres¬ 
sions can easily be found (see Appendix I). As it turns out the result can be analyticaly 
expressed, like the leading corrections, in terms of logarithms, dilogarithms and the two 
functions g{x) and f(x, 1). The function g{x) stems from the on-shell counterterms, while 
f(x, 1 ) from the two loop scalar integrals; they too can be expressed in terms of logarithms 
and dilogarithms. Explicit intermediate results for each of the set of graphs of Fig.2-Fig.7 


are given in the next section. 
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5. Calculations and results 


In this section we present the explicit results of the two loop one particle irreducible 
contributions A^ 2 -* to the universal part A of the p parameter. We neglect the contributions 
of the light fermions and consider only the effect of a single fermion doublet (f, 6), with large 
mass splitting. The mass of the lightest partner in the doublet has been set equal to zero 
(nib = 0) from the beginning of the calculation; this has produced no mass singularities. 
We decompose A^ 2 -* in three parts as : 


A (2) 


A 


( 2 ) 

lead, 


+ A 


( 2 ) 

sub 


+ 


(5.1) 


^ (2) 

where A) e A are the graphs with scalars and fermions that contain leading contributions 

✓s / o \ 

proportional to as well as m 2 , A)^ are the graphs with fermions and gauge bosons that 

✓s / o \ 

contain contributions of order m 2 only, and finally A bos is the pure bosonic contribution 
which is independent of m t and in our approximation is zero. The results are given in 
terms of two variables h = Mfj/m 2 , and r = M 2 /m 2 , where Mjj is the mass of the Higgs 
boson and m t is the pole mass of the heavy quark (top) in the doublet. 

The leading contributions A) e A are given by [38] 


AlL= N c PR(h,r) 


= N c x z t { 23 — 4 h + 


— (1 — r) 2 ln(l — r) + 
2 h 2 


2 h 

h — i 


7T 


— Hr H- r - 2 - 


2 

2 h 2 


Vhg(h) + —\/rg(r) + r 2 In 2 r 


—6 — 2h - 12r + -r 2 + 

2 (h — r) l 


, (1 + h) 

+ 2 h V + 

h — r 


4 — r 


In i 


+ 

+ 

+ 


, h 2 

—Ah -\~ — — /, xo 

2 (h — r ) 2 


- 2 h 


(2 + h ) 

h — r 


In h 


-(1 - h) 2 + 2h 


2 , , 0 ( l -^) 3 + ( l -^) 2 ) 


(h — r) 2 


+ 2 - 


h — r 


Li 2 (l — h) 


-11 + 10 h - 2 h 2 - 2 hr + 16r - 5r 2 - 2 h 


(1 ~hf Jl-hf + il-hf 

(h — r) 2 


- 2 


h — r 


Li 2 (l - r) 
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+ 

+ 


6 — Ah + h 2 + 2h 


2 ( 4 - h) _ 2h l - (3 - h) 
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(h — r ) 2 


h- 


26 - I4h + 2 h 2 + 2 hr - 30r + 7r 2 - 2 h 2 ^ -- 2 2 -^ + 


/(M) 


( A ~ h ) nl. 1 + ( 3 - h f 


(h — r ) 2 


h — r 


4 — r 


f(r, 1) 

(5.2) 


where f(x, 1) is a function stemming from the two-loop scalar integrals [39-40], [ 6 ] 


f(x , 1) = - 


1 


V 4 - y 

= —4 In 2 

2 


y + 2Li 2 (C) + ^ln 2 C 


--Cl 2 ( v ) 


Vv - 4 

and g(x) is a function that originates from the on-shell counterterms 

g(x) = y/4 — x [i t — cf>\ 0 < x < 4 


4 < x 
x = 4 
0 < x <4 


(5.3) 




= 0 

= \/x — 4 In ( —C) 
The variables y, C, and <f> are defined as in [ 6 ] 
4 


x = 4 
4 > x 


(5.4) 


Vi -y -1 , j „ . / 

S = — , - , and d> = 2 arcsm 


V x ’ V 1 - V + 1 r V 2 

By taking the r —s- 0 limit in Eq. (5.2) one can recover the result of Eq.(12) of [ 6 ] where 
the leading contributions were calculated in the approximation Mz = Mw = 0. In the 
analytic formula of Eq. (5.2) the apparent mass singularities cancell in the relevant limits 
h —> r and r —> 4. For example, in the r —> 4 limit the terms that contain r — 4 in their 
denominator give 4(2 In r + /(r, l))/(4 —r) —>■ (4/3)( — 1 + ln 2) . In the limit h —>■ r Eq. (5.2) 
reduces to 

2 

A^Jad = N c x 2 t 119 — 17r -|-r 2 — (2 — r)y/rg(r) + r 2 In 2 r — (1 — r) 2 ln(l — r) 

h=r *• 3 

— 2r(10 — r) In r + 2(— 6 + 14r — 5r 2 )Li 2 (l — r) + 4(11 — 13r + 3r 2 )/(r, 1) j 


Finally taking the limit r —> 0 one obtains the result of [4] 


A 


( 2 ) 

lead, 


= —N c x 2 (27r 2 — 19) 


(5.5) 


(5.6) 


h=r= 0 
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where the leading contributions of a very heavy quark were calculated with all other masses 
neglected. 

^ ( 2 ) • 

We now present our result for AW in full analytic form. We decompose the result as 


A ( 2 ) — N x 2 

^sub — ly c-x t 


aJ 2) + A^ + A^ + A^ + AlA + A(A + A 


( 2 ) 


( 2 ) 


( 2 ) 


.( 2 ) 


.( 2 ) 


.( 2 ) 

6ct 


(5.7) 


and give the analytic form of each intermediate result A) ’ as well. Each A) ’ equals the 
contribution coming from the graphs shown in Fig.2-Fig.7 respectively. This grouping of 
graphs is dictated by their topology, the particles they contain in the loops, and the WI 
of Eq. (2.14). 

The graphs of Fig.2 give 

h 2 , 9 , / 1 h 2 


Aj 2) =4 r‘ 


7r 


12 4 (h - r ) 2 

(1 ~hf 


In 2 h+(-- + 


4 4 (h — r ) 2 


In 2 r 


+ 


(h — r ) 2 
(2 - h){A-h) 


Li 2 (l - h) + ( -1 + yr —) Li 2 (l - r) 


(h — r ) 2 


(5.8) 


2 (h — r) 2 ’ \ 2 r 2 r (h — r ) 2 

In the difference of the self energies the graphs (2.a ) and (2 .a ) cancell for Mw = Mz 
and need not be computed. In the h = r limit the above formula reduces to 


A 


( 2 ) 


= 4r 


h=7 


7r 2 

- 7^4 


1 1 1 2 # , r — 2 /1 

i .--mr- In r —Li 2 (1 —r)H-—I-:-- 

12 4(4- r) 4 v 2 r \2 r(4 - r) 


1 


/(l,r) 

(5.9) 


For the uninteresting value r = 4 the above result as well as the results that follow are all 
regular as can be seen from the explicit values of /(4,1) and / (4,1) given in the Appendix. 
The contribution of the graphs of Fig.3 is given by 


A^ 2) = 2 r 


1 + 


1 


+ 1 - 
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h — r \ h — r 
h( 4 — h) h 2 — Qh + 6 


(h — r ) 2 h — r 
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(5.10) 
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When h = r this reduces to 


A 


( 2 ) 


= 2 r 


h =7 


2 + ( 2 + ^ lnr + (2r-3)Li 2 (l-r)+ (l - 2r - - + ——- ) /(l,r) 


4 — r 


r r(4 — r) 


(5.11) 


The graphs of Fig.4 give 


A< 2> =r 


12 7T^ f 1 3 \ 

--In r + r In 2 r + — r + 4(2 + r)Li 2 (l — r) + 2 I — 21-r — -7 -J f(r, 1) 

(5.12) 


I 


Graphs (4.a) and (A.a) as well as (4.c ) with ( 4 . 7 ) are equal for Mz = Mw and cancell in 
the difference. 

From the graphs of Fig.5 we obtain 


Af } = 2 r 


+ 6 -1 + 


1 


4 — r 


2tt 


In r — 2r In 2 r - r — 6 Li 2 (l — r) + 14 — 4r- 1 -) /(r, 1) 
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(5.13) 


The graphs of Fig .6 give 

— 2r + r 2 (4 — r)- \- 2 
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— 2r In r + r (4 — r) In 
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From the graphs of Fig.7 we obtain 


A< 2) =r 


3 47 7 r 2 

-b 6 £ m + — - 30r - r(4 - 7r) — 

c ^ o 


+ 10 3r + -- J In r — r (4 — 7r) In 2 r+ 

2( —7 + 14r - r 2 )Li 2 (l - r) + 2 27 + —-- - 34r + 8 r 2 ) /(r, 1) 


(5.15) 


where = 7 # + ln(7rm 2 ). 

This last set of graphs requires counterterms which are solely due to the one-loop mass 
renormalization of the top quark. We perform the renormalization on-shell. The fermion 
two-point function Tf is written as 


Tf(p) = i{jt> - m of ) - iZf(p) 


(5.16) 
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where Ej is the fermion self-energy function. The superscript G will denote the gauge 
boson’s contribution to Ej shown in Fig.9. It is given by 


Sf(p)=)SPiS° I (p) (5.17) 

with 

S?i(p) = [Bi(p 2 .m?,M|) + 2B,(p 2 .0,M5,)] (5.18) 

c 

B\ (p 2 , m 2 , M 2 ) = — — Fi (p 2 , m 2 , M 2 ) , (5.19) 

where 

1 1 dz ■'y^ 

S = - - ~ ln7r , and Pr,l =—-— (5.20) 

and 

F n (p 2 , m 2 , M 2 ) = / ln((l — x)m 2 + xM 2 — x(l — x)p 2 ) . (5.21) 

Jo 

The mass counterterm is defined by 

m ot = m t - Sm t (5.22) 

and is determined by the on-shell renormalization conditions 

r*(p>t(p)|y =mt = 0, u t (p)Y t (p) |y =mt = 0 (5.23) 

- = (P =m t ) 

m t 2 

= 32^ “!(l “ [T£ + - 1 + (5.24) 

5r 2 — 12r (1 — r) 2 r — 2 ^ 

-lnr -|-m|l — r|- v r 9y) 

8 2 8 

^ (2) 

This renormalization of the top mass gives rise to two kinds of counterterms for A), j. 
The mass-insertion counterterms shown in Fig.7(ctl,ct3) and the vertex counterterms 
shown in Fig.7(ct2,ct4). The mass insertion counterterms diverge as 1/e, which cancells 
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in the difference, and they give a finite contribution to A), j. In terms of scalar integrals 
this contribution equals 
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6ct(mi) 
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where henceforth m = m t . The vertex counterterms turn out to be just the one loop 
fermion graphs multiplied by the factor qjp. In terms of scalar integrals they are given by 

,( 2 ~ D ) M 
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6 ct(v) 




. g 2 m 2 \ 8m 


(27t) 4 M 2 [ m \ [ 
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+ (mm) 
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To see how this comes about, we note that upon renormalization the (fiifiif) (and similarly 
the Xtjitf) ) vertex will get modified to 


M o 


+ , /2 

9y (M + SM) * 


where we have ommited wave function renormalization of the fermion fields, since it will 
cancell agianst the corresponding wave function renormalization of the fermion propagators 
inside the loop. Then, at one loop the modification reads 


8m 8M 1 
m M 


1 + SZg + — — + T,^ Z <t> ~ 1 + 


8m 

m 


where in the last step we have used Eq. (2.20). Finally the result for the counterterms is 


A (2) -r 

^6 ct 


3 5 / 5 

- 6£ m H-3r — r ( 6 - r)lnr — 2(1 — r) 2 ln(l — r) 

e 2 V 2 


1 


+ r( r “ '2)Vrg(r) 


(5.27) 


We note that after the inclusion of the counterterms, the terms proportional to 1/e 
( 2 ) 

present in Ag J cancell, and the final result emerges finite as expected. In the limit where 
the mass spliting in the doublet is zero, which in our case means m t = mj = 0 , all of the 
above expressions vanish, as can be seen from the asymptotic expressions of the functions 
g(x), f(x, 0 ), and f(x, 1 ), if we take x —> oo (for x = r or h). 
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Adding together Eq. (5.8) to Eq. (5.27) we report as our final result 
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In Figure 10 we show together the functions R(h } r ) and S(h } r) } which describe re¬ 
spectively the leading and subleading contributions, for M/m t = 0.5 or r = 0.25 and a wide 
range of values for the mass of the Higgs boson, Mjj, where h = Mfj/m 2 . The function 
R(h } r ) of the leading contributions grows asymptoticaly in the way described in [ 6 ]. On 
the other hand, the function S(h } r ) of the subleading contributions becomes independent 
Mjj for a very heavy Higgs boson (as pointed out in [7]). For the phenomenologically 
interesting range of values for Mjj , the two functions have opposite sign. In Table 1 we give 
numerical results for A^ 2 -* for different values of the ratios M/mt and Mu/mt in units of 
N c x 2 . The first entry of each column corresponds to the leading contribution while 

the second gives the total correction ^) ead + We notice that for a light Higgs boson, 

where A^ 2 -* is small, the two contributions ^\ ead and A) j are comparable in magnitude, 
as can be seen from the entries in the first three lines of the first four columns. We also 
observe that for a light top ( M/m t = 0 . 6 , fifth column) the two corrections come with the 
same negative sign for Mjj A mt- On the other hand, for the largest part of the parameter 
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^ ( 2 ) 

space of rrit and Mjj, the subleading contributions are approximately 22-27% of the 
^( 2 ) 

leading 7S. lead contributions, which is what one would naively expect. Finally, in Figure 
11 we give the two loop correction A^ 2 -* as a function of Mu for M = 91.19 GeV and 
rrit = 175 GeV. 


6. Conclussions 

In this paper we showed that in the framework of the PT one can define a universal 
part of the p parameter, which satisfies all necessary field theoretical requirements, for both 
leading and sub-leading two-loop corrections. Most noticeably, the PT universal part is 
by construction independent of the gauge-fixing parameter, and, at the same time, process 
independent; furthermore, by virtue of the PT Ward identities it is also UV finite. 

We have calculated the two loop contributions at the subleading order (9(G 2 m 2 M^), 
in the limit Mw = Mz, and s = 0. Their relative size was found to be around 25% with 
respect to the leading ones. From the technical point of view, the computation involved 
self-energy garphs only, which exist in closed analytic form. The computational part was 
significantly fascilitated by the PT Ward identities, relating the self-energies of the gauge 
bosons (WW and ZZ) to the corresponding self-energies of the would-be Goldstone bosons 
(</></> and XX). These Ward identities, which are valid for both leading and sub-leading 
contributions, reduce the task into calculating Goldstone boson self-energies only. 

Having laid out the framework of how such calculations should proceed, it is straight¬ 
forward to compute the two loop corrections to the universal part of p defined via the PT, 
for the full SU{2)l X U{1)y standard model. Results in this direction will be presented 
elsewhere. 
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8. Figure captions 

Figure 1 : The £ q dependent part of A in the background field gauges at one loop 
(for Mh = 300 GeV and in units of g 2 / (IQtt 2 )) . 

Figure 2 : Subleading two loop graphs containing self energy insertions of a Higgs 
boson or a gauge boson. 

Figure 3 : The subleading fermionic triangle graphs that contain a Higgs boson. 

Figure 4 : The rest of the subleading graphs containing bosonic self energy insertions. 

Figure 5 : The rest of the subleading triangle graphs. 

Figure 6 : Subleading fermionic bubble graphs with a vertex correction. 

Figure 7 : Subleading fermionic bubble graphs with self energy insertions and their 
relevant counterterms. 

Figure 8 : Subleading graphs that contain a mixed self energy insertion. The graphs 
that correspond to the scalar self energy </></> vanish for sin 6w = 0. 

Figure 9 : The gauge boson contributions to the one loop top self energy. 

Figure 10 : The functions R(h } r) and S(h } r) describing respectively the leading and 
the subleading contributions to A, for r = 0.25. 

Figure 11 : The universal two loop A^ 2 -* correction to the p parameter as a function 
of the Higgs boson mass, for different approximations, in units of N c x 2 . Dotted line : the 
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leading result in the approximation Mw = Mz = 0, [5], [6]. Dashed line : The leading 
contribution in the approximation Mw = Mz = M [7]. Solid line : the total correction, 
leading and subleading, for Mw = Mz = M. ( m t = 175 GeV, M = 91.19 GeV). 

Table 1 : Numerical values for A^ 2 ) in units of lV c x 2 , for different values of \fh = 
Mu/mt and yjr = M/m t . The first entry of each column corresponds to the leading 
contribution A lead in the approximation Mw = Mz = M, while the second gives the total 
result AG) = Aj ead + A { s J b in the same approximation. 
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Appendix I. 

Scalar two loop integrals 


The scalar two loop integrals at zero external momentum are defined as in [39] and 


[40]: 


. . . , . . d D kd D £ . . 

(m 1 ...m 1 \m 2 ...m 2 \m 3 ...m 3 ) = / — - - 2 - (1. 1 ) 

11 [k 2 — m\] [£ 2 — m\] [(k + £) 2 — m\] 


V\ V 2 ^3 


Obviously 


(mi...mi | m2...m2 \ m3...m3) = (m2...m2 \ mi...mi \ m3...m3) = (m3...m3 \ m2...m2 \ mi...mi) 


V\ 


V2 


V 3 


V2 


V\ 


V 3 


V\ 


V2 


V\ 


Note that, within each of the three groups of masses in Eq. (1.1), all masses are equal. Any 
other integral, for which different masses appear within the same group, can be brought 
in the form of Eq. (1.1) by use of partial fractions 


1 1 _ 1 1 1 

k 2 — m 2 k 2 — m\ m 2 — m\ k 2 — m 2 k 2 — m\ 

All of these integrals can be calculated from the knowledge of a single master integral, either 
by use of recurrence relations, derived through integration by parts, or by differentiation 
with respect to one of the masses. For example 


( 1 . 2 ) 


(MMMIM, |M 2 ) = |M 2 ) 

and 

(MM\MiMi\M 2 ) = -^(MM\Mi\M 2 ) . 

For our purposes, we have found the following two recursion relations very useful 
(M 1 XM 3 M 3 ) - (M 3 )(MiMi) + (M 2 ) ((MiMi) - (M 3 M 3 )) = 

(Ml - Ml - M 3 2 )(M 3 M 3 |Mi|M 2 ) + (Ml + M| - Mf)(MiMi |M 2 |M 3 ) 


(/.3) 

(1.4) 

(1.5) 
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( 1 . 6 ) 


(D — 3)(Mi |M 2 |M 3 ) =M 1 2 (M 1 M 1 |M 2 |M 3 ) + M|(M 2 M 2 |M 1 |M 3 ) 

+ M 3 2 (M 3 M 3 |M 1 |M 2 ) , 

The integrals of the type (MM), are one loop integrals and are defined at the end of this 

Appendix. To prove Eq. (1.5) we consider 

2k ■ q 


Q( n 2 '\ = / r] D krl D f _ 

J [(£ + ?) 2 — m i] (k 2 ~ m 2 ) [(k + q + k ) 2 — m 2 ] 

which is zero as can be easily seen by the shift of variables £ —> £ + q. Then from 


( 1 . 7 ) 


dS(q 2 ) 


dq 2 


,,4M „=° 

q z = 0 \ / g=0 


( 1 . 8 ) 


one obtains Eq. (1.5). To prove the identity of Eq. (1.6), which has also been given in [39], 
we use 


L(q 2 ) = / d D kd D £ 


d 

dkv 


k. 


= 0 


(/.9) 


( k 2 - Ml)(P - Ml)[(k + Cf- - M|] 

and Eq. (1.5). 

In [39], the master integral (MM|M 3 |M 2 ) has been calculated . In [40] the same 
integrals have been discussed in a different mathematical framework and the master integral 
(M|Mi |M 2 ) is given. In what follows we use the notation of [39]. With D = 4 — 2e the 
result for the master integral is : 


(MM|Mi|M 2 ) = 7r 4 

where 


'^2 “ 2e^ ~ ~ k 2 M - - - — - /(a, b) 


(I- io) 


M 2 M 2 

k M = + In 7 tM 2 , a=^ 6=^. 

The real function /(a, b) is symmetric in its arguments and defined as 


d 


f(a,b) = j dx Li 2 (l — ji 2 ) — i r 2 ln/i 2 


1 — ji 2 


d = 


ax + 6(1 — x) 
t(l — x) 


with Li 2 the dilogarithm function 


U 2 (x) = - [ X dy Hl V) 

Jo y 


(i- ii) 
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The explicit form of /(a, b ) is 


1 


f(a.b) =-In a In 6 + , 

; 2 yf L 6 


1 — a — 6 r 7T 2 _ . / x 2 \ , T . ( y 2 


— + Li 2 - — +Li 2 

yi / V *1 


+ I ln ^ + i h ^ + i ln *£! 

4 yi 4 Xl 4 


yi 


i l n 2 ^ 

4 


y2 j 


(/. 12 ) 


with 


V" = \Jl — 2 (a + 6 ) + (a — 6 ) 2 , xi ;2 = - 1 + 6 —a 


’ yi ’ 2 2 L 


1 + a — b ± -s/ - 


It turns out that all of our results can be expressed in terms of /(a, 0 ) and /(a, 1 ). 
These functions are explicitly given by : 


7r 


/(a, 0 ) = LWl — a) = -LWa) — In a ln(l — a) , 

6 


(7-13) 


and 


f(a, 1 ) = - 


1 


V 4 - y 

= —4 In 2 

2 


Vy - 4 


-y + 2 Li 2 (C) + ^ln 2 C 


Cl 2 (p) 


where Cl 2 (x) is the Clausen function 


4 < a 
a = 4 
a > 4 


( 7 - 14 | 


Cl 2 (:c) = Im[Li 2 (e")] , 


( 7 . 15 ) 


and the variables y, C, and <f> are defined as in [ 6 ] 
4 


y = 


y V 1 - y - 1 ! i o ■ ( V x 

S = —-- , and <p = 2 arcsm - 


•s/i — y + 1 

The derivatives of /(a, 1) and /(a, 0) are given by 


df(a , 1 ) 

da 


4 — a 


In a H-/(a, 1 ) 


= --(l + 21n2) , 


a / 4 , 
a = 4 


df(a, 0 ) In a 


da 


1 — a 
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We will also need the following derivative 


a/(a, b) 


db 


6=1 


4 — a 


In a + ( 1-) /(a, 1) 

a ' 


1 4 , 

=-In 2 , 

3 3 


, « / 4 

a = 4 . 


which we obtain directly from Eq. (1.12). 

The following relationships are useful in expressing the results in terms of only these 
two functions 


/(-, 0) = -/(a,0) - ^ln 2 a 

a 2 

/(-,-)= - fl --)/(«, 1) - ^ln 2 a 

a a V a / 2 


(7.16) 


In our calculations we have encountered two loop integrals with three, four, and five 
propagators. All integrals with three propagators are converted to integrals with four 
propagators using Eq. (1.6). Most of the integrals with four propagators can be directly 
obtained from Eq. (1.10). For the integrals (00|m|M) and (00|0|M) we use Eq. (1.5), and 
write them as 

(1 — r)(00\m\M) = —— (1 + r)(MM\m\0) 




and 


with r = M 2 /m 2 as usual. 


(00|0|M) = —(MM|0|0) 


The integrals with five propagators that we encountered are the following 

1 d 

(MMM\m\m) = — ( MM\m\m ) 


7r 


2 M 2 


1 21nr 4 .. . 

- + 1 - 2£ m + - b —- -/(r, 1) 

e 4 — r r( 4 — r) 


(MMM\m\0) = I^A_(MM|m|0) 


7r 


2 M 2 


1 „ In r 

b 1 — 2 £m + z 
e 1 — r 


(7.17) 


(MA7M|0|0) = 


7T 


2M 2 


—b 1 — 
e 
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1 d 
= — —— 

2 dm 2 


7r 


2 m 2 

r 4 


fi=m 

1 „ df(a.r), 

e da l‘=i 


_l_ r 


dj(l,r) 

dr 


7r 




2 m 2 

d 

dfi 
4 2 


- + 1 - 2C - ^ Inr + —-/(r, 1) 

e 4 — r r(4 — rj 


7T 4 d/(r,6) 


fi=m 




db 


(I- 18 ) 


6=1 


7T 


m 2 4 — r 


In r + ( 1-] /(r, 1) 


d d 7f 4 df(v ll 

(MM|mm|0) = —— -(MM\m\0) = = —- ’ 


dm 2 
7T 4 In r 
m 2 1 — r 




dr 


Besides the two loop integrals our results also contain products of one loop integrals. 
These integrals can all be computed from the integral (m) through differentiation. 




i 


1 


d r ‘ 


[k 2 — m 2 ] n (n — 1)! d(m 2 ) 

n 

The expression for (m) up to (9(e) is 
(m) 




i 


k 2 — m 2 


= Z7T m 


+ l-k m + e(l + — - £ m + ~k 2 m 


(7.19) 
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Appendix II. 

Feynman rules in the background field gauges 


In this appendix we give the Feynman rules that are relevant for the calculation of the 
subleading two loop contributions to A. The rules are given for the full standard model, 
and the calculations in the text were carried out in the approximation : Mw = Mz } s = 0, 
mj = 0. 
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